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SPATIAL  REALIZATIONS  OF  SHIP  RADIATION  AND  DIFFRACTION 
WAVES  IN  RANDOM  AMBIENT  SEA  WAVES 


1.  INTRODUCTION 

The  calculation  procedure  as  well  as  numerical  examples  for  the  radiation  and  diffraction  wave 
patterns  caused  by  a  ship  at  zero  forward  speed  excited  by  a  single  sinusoidal  ocean  wave  (of  given 
frequency  and  heading  relative  to  the  ship)  have  been  given  by  the  present  author  in  a  previous  report 
[1].  Radiation  waves  refer  to  those  caused  by  the  six  oscillatory  modes  of  motion  of  the  ship,  while 
diffraction  waves  refer  to  the  scattering  of  the  incident  wave  by  the  ship  hull,  taken  to  be  stationary. 
It  was  pointed  out  in  [1]  that  most  of  the  calculations  involving  a  ship  excited  by  ambient  waves  have 
been  largely  concerned  with  near  field  flow  quantities  such  as  motions,  forces,  and  wave  profiles. 
However,  with  the  introduction  of  modem  day  observation  techniques  such  as  airborne  and  space- 
borne  Synthetic  Aperture  Radar  (SAR)  systems  and  other  methods  of  remote  observation  such  as 
infrared  (IR)  and  microwave  (MR)  radiometry,  the  interest  has  been  enlarged  to  include  the  far  field 
wave  pattern.  A  number  of  recent  studies,  such  as  [2,3],  have  focused  on  the  far  field  features  of  the 
Kelvin  wake,  the  wave  system  behind  a  nonoscillating  ship  advancing  at  constant  forward  speed. 

The  present  work  was  undertaken  to  obtain  estimates  of  the  far  field  radiation  and  diffraction 
waves  of  typical  ship  hulls  excited  by  random  seas.  Essentially,  the  previously  developed  calculation 
procedure  for  a  single  sinusoidal  wave  is  generalized  to  a  random  sea,  composed  of  a  series  of  com¬ 
ponent  waves  of  different  frequencies  and  headings. 

The  report  starts  with  a  summary  of  the  calculation  procedure  for  a  single  wave  component  [1]. 
It  is  shown  that  the  wavemaking  due  to  all  the  sources  on  the  hull  surface  is  expressed  in  terms  of  the 
concise  and  convenient  Kochin  function.  In  this  summary,  as  well  as  throughout  the  report,  the  ship 
waves  and  ocean  waves  are  conveniently  given  in  compact  complex  notation.  This  makes  it  unneces¬ 
sary,  for  example,  to  explicitly  give  the  phase  of  the  ship  waves  with  respect  to  the  ocean  waves. 
The  generalization  of  the  single  wave  procedure  to  the  multiple-wave  random  sea  is  discussed  in 
detail.  This  includes  the  technique  for  selecting  the  frequencies  and  headings  for  modeling  the  ran¬ 
dom  sea  itself,  the  superposition  of  the  ship  waves  due  to  each  sea  wave  component,  and  the  calcula¬ 
tion  of  the  ship  and/or  sea  waves  at  any  arbitrary  point  on  the  ocean  surface.  The  formulation  is 
given  for  the  slopes  as  well  as  the  elevations  of  these  waves.  The  procedures  required  for  adapting 
the  transfer  function  and  Kochin  function  output  at  the  selected  frequencies  and  headings  of  the  Ship 
Motion  Program  SMP  [4]  to  those  modeling  the  random  sea  are  discussed  in  some  detail. 

The  preceding  results  are  also  generalized  to  the  case  of  a  twin  hull,  by  essentially  reflecting  the 
Kochin  functions  for  a  single  hull  about  the  vertical  plane  of  symmetry  of  the  twin  hull.  The  errors 
incurred  by  using  the  results  of  SMP,  which  is  for  a  monohull,  to  generate  waves  for  a  twin  hull  are 
indicated. 

Spatial  realizations  are  then  presented  for  the  ship  and/or  sea  wave  elevations  and  slopes  for 
four  ship  hulls  in  random  seas  corresponding  to  Sea  States  2,  4,  and  8.  The  spatial  realizations  are 
presented  as  three-dimensional  plots,  contour  plots,  and  two-dimensional  plots  of  wave  elevations 
along  a  given  line.  Spatial  realizations  such  as  these  do  not  appear  to -have  been  previously  given. 
The  report  concludes  with  a  summary  of  the  principal  findings. 

Manuscript  approved  December  15,  1986. 
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2.  FORMULATION  FOR  SINGLE  EXCITING  WAVE 
2.1  Monoksll  Case 


Adopting  the  usual  assumption  of  linearized  potential  flow,  the  problem  essentially  consists  of 
two  parts:  the  evaluation  of  the  potential  flow  on  and  near  the  ship  hull  and  the  extrapolation  of  this 
near  field  solution  to  obtain  the  far  field  waves. 

The  near  field  problem  consists  of  finding  a  potential  function  <t>  which  satisfies  Laplace’s  equa¬ 
tion,  the  free  surface  condition,  the  radiation  condition  of  outward  progressive  waves  at  infinity,  and 
the  kinematic  condition  of  no  flow  through  the  ship  hull.  In  the  commonly  used  Green’s  function 
approach,  the  procedure  is  to  place  a  series  of  singularities  on  the  hull  surface  which  satisfy  the  first 
three  conditions,  and  then  to  determine  their  strengths  from  the  kinematic  condition  on  the  hull  sur¬ 
face.  In  the  present  work,  these  strengths  are  calculated  by  using  the  Ship  Motions  Program  SMP 
[4],  which  neglects  longitudinal  interactions  between  ship  cross  sections,  thus  reducing  the  three- 
dimensional  problem  to  a  series  of  simpler  two-dimensional  problems. 

The  coordinate  system  used  in  the  calculations  in  SMP  (and  also  used  in  the  present  work)  is 
shown  in  Fig.  1.  The  origin  lies  on  the  undisturbed  free  surface,  at  the  longitudinal  center  of  gravity, 
and  in  the  vertical  plane  of  symmetry.  The  x-axis  is  directed  forward,  the  y-axis  to  port,  and  the  z- 
axis  upwards.  In  this  coordinate  system,  following  and  head  seas  correspond  respectively  to  headings 
of  0  and  180  degrees.  It  should  be  noted  that  this  coordinate  system  and  sea  heading  convention  do 
not  correspond  to  those  used  for  input  and  output  in  SMP. 

To  first  order,  the  strengths  of  the  sources  obtained  for  the  near  field  may  also  be  used  to  obtain 
the  far  field  waves  [5].  However,  the  sources  themselves  must  be  changed  from  those  for  two  dimen¬ 
sions,  Gw,  to  those  for  three  dimensions,  G3D.  In  the  present  work,  the  resultant  wave-making 
capability  of  the  sources  distributed  over  the  hull  surface  is  obtained  in  terms  of  the  Kochin  function, 
which  represents  a  weighted  integral  of  the  source  strengths  over  the  hull  surface. 


Wehausen  and  Laitone  [6]  give  the  following  asymptotic  expression  for  the  far  field  wave  eleva¬ 
tion  fy  of  the  radiation  wave  due  to  an  oscillation  ^  of  unit  amplitude 


(1) 


where  j  —  1,. . 6  denote  the  six  modes  of  ship  oscillation 

k  —  —  —  —■  is  the  wave  number 
8  X 

03  is  the  circular  frequency  of  the  exciting  wave 
X  is  the  wavelength  of  the  exciting  wave 

R  =  'lx2  +  y2  denotes  the  radius  in  the  horizontal  plane,  measured  from  the  origin  shown  in  Fig.  1 
0  =  tan”'(y/x)  is  the  direction  in  the  horizontal  plane,  measured  from  the  x  -axis 
Hj  is  the  complex  conjugate  of  the  Kochin  function  Hj 
GjD  is  the  asymptotic  expression  for  the  three-dimensional  source  potential  at  large  R 
The  Kochin  function  Hj  is  given  by 

Hj  ( k , 0 )  -  -  JJ  +  '*,*cos<’  +  '>*inB)'dS  (2) 

s 
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where  yj  is  the  hull  surface  distribution  of  singularity  strengths  calculated  by  SMP  for  the  jth  mode 
of  ship  oscillation;  are  respectively  the  x,y,z  coordinates  of  the  source  located  on  the  huli  sur¬ 
face.  It  is  shown  in  [7,8]  that  at  distances  R  greater  than  0.3  X,  there  is  relatively  little  difference  in 
die  wave  elevations  calculated  by  the  asympototic  and  complete  formulas.  In  view  of  this,  Eq.  (1)  is 
used  to  calculate  tj  for  R  5:  0.3  X.  For  R  <  0.3  X,  the  coefficient  ^2k/v R  is  simply  set  equal  to 
^2k/ir(0.3\y,  otherwise,  —  oo  as  R  —  0. 

The  diffraction  wave  fR7  for  an  incident  wave  of  unit  amplitude  is  given  by  the  sum  of  the  ver¬ 
tical  and  lateral  components  £R7y  and  %R7L 

?«7  “  $R1V  +  ftf7I 

=  Re  -  \f~^  (H7V  (k  ,6,(3)  +  H7L  (k  ,0,(3))  em  +  -  ~  r/4)  (3) 

g  2vK 

where  £  is  the  wave  heading,  H7V  and.  H7L  are  respectively  the  complex  conjugates  of  the  symmetric 
(vertical  plane)  and  antisymmetric  (lateral  plane)  components  of  the  diffraction  Kochin  functions 
H7V  and  H7L ,  which  are  evaluated  from  Eq.  (2)  with  the  singularity  strengths 

■ftv  -  -7* -*«•“»  (4a) 


77t  =  i  sinj372«-<*£c0!6 


For  the  case  where  the  ship  oscillation  is  not  of  unit  amplitude  but  is  the  response  to  a  wave  of 
amplitude  a  arriving  with  phase  <fs  relative  to  the  origin  of  the  coordinate  system,  in  Eq.  (1)  must 
be  multiplied  by  ae **  and  the  transfer  function  7}(u,/3),  while  must  be  multiplied  by  ae,4>  .  In 
the  present  work,  for  the  sake  of  computer  programming  convenience,  the  resultant  of  the  radiation 
and  diffraction  waves,  ^ ,  is  calculated  by  one  single  expression,  as  follows 

SrdM,®  =  Re  E  cte^TjM)^  +  a e'+ifay  +  ?R7L) 

8  8 

=  Re  ete'+Z  Tjfafitj  =  Re  £  f RDj  (5) 

j= i  J  L>-» 

where  £7  =  Zriv 


fg  “  Zril 

Tj(o),(3)  a 


j  =7,8 


2o2  Twin  Hull  Case 


In  order  to  obtain  an  estimate  of  the  wave  patterns  due  to  twin  hull  ships,  the  above  formulation 
for  monohulls  has  been  generalized  to  approximately  model  the  twin  hull  case.  The  procedure  con¬ 
sists  of  summing  the  preceding  expressions  for  the  Kochin  function,  which  are  formulated  for  a 
monohull  with  vertical  centerplane  at  y  =0,  for  two  hulls  with  identical  singularity  distributions  at 
y  sas  y(,  and  y  =  Y2.  From  Eq.  (2),  the  resultant  Kochin  function  HT)  for  the  two  hulls  is  given  by 
the  following  sum  for  two  monohulls 


UTj  =  Hjiy  =  Y{)  +  Hj(y  =  Y2) 

,,  ,  ikY',  sin 5  ,  /fcK,sin0, 

-Hj{e  +  e  ■  ), 


j  -  1,2,... 6, 7, 8 
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For  the  particular  case  where  the  two  hulls  are  symmetrically  located  2  Y  apart  about  y  -  0  Hr  is 

givers  by  ’  - 

HTj  -  Hj  (e-ikY™6  +  eikY™*)  =  Hj  2  cos  (kYsind)  (7) 

The  above  procedure  gives  the  wavemaking  of  two  hulls  which  do  not  affect  each  other’s  motion 
and  whose  waves  completely  pass  through  the  adjacent  hull.  In  the  actual  twin  hull  case,  there  is 
both  hydrodynamic  and  dynamic  coupling  between  the  motions  of  the  two  hulls  and  the  channel 
between  the  two  hulls  is  an  area  of  mutual  wave  reflection. 

In  the  case  of  roll  motion,  Jones  [9]  experimentally  shows  that  the  motions  decrease  with 
increased  hull  spacing,  largely  due  to  the  increase  in  roll  moment  of  inertia.  For  a  three-fold  increase 
in  hull  spacing,  Jones  obtains  a  two-fold  decrease  in  the  roll  motion  transfer  function  at  the  resonance 
frequency,  with  less  difference  at  other  frequencies.  It  is  of  interest  to  note  that  for  a  ship  with  zero 
forward  speed,  the  case  considered  in  the  present  work,  there  is  little  change  with  hull  spacing  in  the 
vertical  plane  pitch  and  heave  transfer  functions.  Thus,  the  present  approach  overestimates  to  some 
extent  the  roll  motions  but  may  be  reasonably  accurate  for  the  vertical  plane  motions. 

Lee  [10]  accounts  for  the  hydrodynamic  interaction  between  hulls,  within  the  strip  theory 
approach,  by  modifying  the  usual  two-dimensional  Green’s  function  to  consider  the  vertical  center- 
plane  between  the  two  hulls  as  a  plane  of  symmetry.  This  leads  to  a  region  of  completely  trapped 
waves  between  the  two  hulls.  It  is  pointed  out  in  [11,12]  that  this  may  lead  to  an  exaggeration  of  the 
interaction  effect  since  some  of  the  waves  inevitably  leak  out  of  the  channel.  It  is  pointed  out  in  [11] 
that  motions  obtained  by  neglecting  the  hydrodynamic  interaction,  as  is  done  in  the  present  work,  may 
lead  to  more  accurate  results. 

Perhaps  the  most  serious  error  in  the  wave  patterns  predicted  by  the  present  approach  is  in  the 
channel  area  of  trapped  waves  between  the  two  hulls.  It  may  be  noted,  however,  that  this  is  a  near 
field  region,  which  is  not  accurately  modeled  by  the  present  far  field  approach. 


3.  WAVE  ELEVATIONS  AND  SLOPES  FOR  RANDOM  SEA  CASE 
3.1  Generation  of  Random  Sea  Elevations 


Following  the  procedure  outlined  in  [13,14],  the  random  sea  elevations  fSE(x,y,t)  are  approxi¬ 
mated  by  a  double  sum  over  wave  number  km  and  wave  heading  n„  of  sinusoidal  components,  as  fol¬ 
lows 


tsE  - 


M  N 

E  E 

m  =1  n  =  1 


O’ 


«'$«.  »'[*„  (*  cos  ft,  +  y  sinpj  -w,/) 


Re 


M  N 

E  E 

n  =  1  n  =  1 


?SEn 


(8) 


where  anm  -  V25(o3m)Aw  ^G(nn  )Ky.  =  aman  (9a) 

S(oim  )  is  the  value  of  the  frequency  spectrum  at  wave  frequency  u)m 
G(ixn )  is  the  value  of  the  spreading  function  at  wave  heading  /t„, 
taken  relative  to  the  wind  direction  &w 

Pn  ~  &n  ~  (9b) 


0n  is  the  wave  heading  measured  with  respect  to  the  x  axis 
and  An  are  increments  in  w  and  n 
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In  the  present  approach,  the  following  three  widely  used  frequency  spectra  are  considered:  Pierson- 
Moskowitz,  JONSWAP,  and  Neumann.  The  spreading  function  G(n)  is  taken  to  be  the  commonly 
used  cosine  squared  function 

G(n)  =  —  cos2  n  —  £  fi  ^  (10) 

7T  2  2 

It  is  shown  in  [13]  that  as  few  as  five  values  of  osm,  chosen  so  as  to  give  equal  values  of  am, 
are  sufficient  to  give  wave  amplitudes  which  approximate  those  of  a  random  sea.  In  [14]  it  is  sug¬ 
gested  that  10  or  more  heading  components  are  needed  to  give  a  visually  random  sea. 

3.2  Ship  Radiation  and  Diffraction  Wave  Elevations  in  a  Random  Sea 

The  resultant  ship  radiation  and  diffraction  wave  f°r  the  random  sea  case  is  obtained  by 
summing  the  ship  waves  for  a  single  incident  wave,  given  in  Eq.  (5),  for  the  above  M  xN 
sinusoidal  components,  resulting  in  the  following  triple  sum 

M  N  8  M  N  8 

ZsH  ~  R®  £  £  amne  =  Re  £  I]  L  f SHmnj  ___  (^) 

j  =  1  m  =  ln  =  ly=l 

The  total  of  the  ship  and  sea  waves,  $>,  is  then  simply  given  by 
fr  =  £se  +  fsw 

M  N  f  8 

=  £  E  fsEmn  +  £  $  SHmnj 

m  =  1  n  =  1  1  j  =  t 


3.3  Wave  Slopes 

The  preceding  formulation  for  wave  elevations  may  be  conveniently  used  to  obtain  expressions 
for  wave  slopes.  Taking  x  and  y  derivatives  of  the  sea  wave  elevations  given  by  Eq.  (8)  results  in 
the  following  formulas  for  the  x  and  y  slopes,  respectively  tSEx  and  fSEy 

M  N 

i*S£x  ~  Re  2  £  ^OS  jln  I" SEmn  (13) 

m  =1  n  «1 
M  N 

%SEy  =  Re  £  £  ikm  sin  Hn  f SEmn  <14> 

m  =  l«  =! 


Expressions  for  the  jc  and  y  slopes  of  the  ship  waves  are  obtained  by  first  taking  the  x  and  y 
derivatives  of  G°D  given  in  Eq.  (1).  Keeping  only  the  lowest  order  term  for  large  values  of  R,  and 
further  noting  the  definitions  given  for  R  and  8  in  Eq.  (1),  dG°D/dx  and  dGyD/dy  are  given  by 


=  ik  cos  9  GjD 
—  ik  sin  9  G°p 


(15) 

(16) 


From  the  summation  given  in  Eq.  (11),  it  follows  that  the  x  and  y  slopes  of  the  resultant  ship  radia¬ 
tion  and  diffraction  waves,  ?SfLx  and  fSHy ,  are  given  by  the  following  triple  sums 

M  N  8 

.  i‘sHx  =  Re  £  £  £  ikm  cos  9  $SHmnj  (17) 

m  =  I  n  =  I  j  ~  I 
M  N  8 

fs'Wv  ”  Re  £  £  £  lkm  Sm  9  f SHmnj  (18) 

m  =  I  n  =  I  j  =  1 
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The  total  of  the  ship  and  sea  wave  .t  and  y  slopes,  £Tx  and  j*7y ,  are  given  by  the  following  sum¬ 
mations 


fs£c  + 

tsHx 

M 

N 

r 

8 

Re  E 

E  i 

k 

COS  fln  r SEmn 

+  km  COS  6  ^  f SHmnj 

(19) 

m  »=1 

n = i 

Vw 

J-\ 

tsEy  + 

?shy 

M 

N 

r 

8 

Re  £ 

E  * 

Ka 

Sin  fln  r SEmr. 

+  km  sin  B  £  tSHmnj 

(20) 

T71  1 

n  =1 

/-»  J 

4.  COMPUTER  IMPLEMENTATION  OF  METHOD 
4J.  Computer  Program  SERDWAV 

The  formulation  has  been  implemented  in  the  form  of  computer  program  SERDWAV.  The  pro¬ 
gram  accepts  from  one  input  file  data  on  the  speed  and  direction  of  the  wind,  the  frequency  spectrum 
to  be  used,  the  number  of  frequency  components,  and  the  types  of  waves  to  be  calculated  (ship,  sea, 
or  total  of  ship  and  sea;  elevation,  x  slope,  or  y  slope).  In  this  same  file  are  specified  data  for  the 
rectangular  grid  over  which  the  waves  are  calculated:  its  origin,  and  the  number  and  spacing  of  the 
points  along  the  ( x',y ')  axes,  which  are  inclined  at  angle  a  to  the  (;t,y)  system  (see  Fig.  2).  In  the 
case  of  a  twin  hull,  this  file  also  gives  the  lateral  offsets  Y{  and  Y2  of  the  vertical  centerplanes  of  the 
two  hulls.  As  a  second  input  file,  the  program  accepts  the  Kochin  functions  Hj(u,d,(3 )  and  transfer 
functions  Tj(oi,0)  calculated  by  SMP  [4]. 

The  program  then  performs  the  calculations  for  the  wave  elevations  and  slopes  given  by  Eqs. 
(g)-(20).  In  the  case  of  a  twin  hull,  the  Kochin  functions  Hj  are  modified  according  to  Eq.  (6).  The 
calculated  waves  are  saved  as  an  output  plot  file  from  which  three-dimensional,  two-dimensional  or 
contour  plots  may  be  generated.  Examples  of  these  plots  are  given  in  the  following  chapter. 

4=2  Adaptation  of  Formulation  to  SMP  Results 

Program  SMP  calculates  the  functions  Hj  and  7}  for  30  discrete  values  of  w  between  0.2  and  2.4 
rad/sec,  and  13  equally  spaced  values  of  6  and  0  between  0  and  180  degrees.  Thus,  it  is  of  interest  to 
describe  the  procedure  used  to  model  the  general  case  of  values  of  «,  0,  and  6  which  may  not  coincide 
with,  or  fall  outside  the  range  of,  the  SMP  results. 

4.2.1  Procedure  for  wave  frequencies  cj 

In  order  to  more  accurately  model  the  lower  sea  states,  where  the  higher  co’s  are  of  interest,  the 
two  highest  frequencies  calculated  by  SMP  were  changed  from  2.2  and  2.4  to  2.4  and  2.8  rad/sec, 
respectively.  This  essentially  consisted  of  changing  one  DATA  statement  in  SMP.  The  actual  values 
of  w  (and  the  corresponding  values  of  wavenumber  k )  appearing  explicitly  in  Eqs.  (8)-(20)  are 
retained.  On  the  other  hand,  the  values  of  Hj  and  7}  are  approximated  by  those  corresponding  to  the 
SMP  values  of  w  which  come  closest  to  the  actual  values. 

4.2.2  Procedure  for  wave  headings  n  and  0 

Since  the  results  of  [14]  indicate  that  a  minimum  of  approximately  10  wave  headings  are 
required  to  give  visually  random  seas,  it  was  decided  to  simply  adopt  the  13  equally  spaced  headings 
between  0  and  180  degrees,  as  used  by  SMP.  Noting  that  the  energy  of  the  wind  is  spread  90 
degrees  on  either  side  of  the  wind  direction  (see  Eq.  (10)),  the  SMP  formulation  is  most  conveniently 
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used  when  the  wind  direction  is  SO  degrees,  i.e.,  in  this  case  the  wave  headings  of  interest  are  simply 
the  0  to  180  degree  cases  considered  by  SMP.  For  other  wind  directions,  however,  Fig.  3  shows  that 
negative  values  of  fi  are  involved.  In  these  cases,  use  is  made  of  the  fact  that  for  the  vertical  modes 
(j  =  1,3, 5, 7)  7}(a>,-/3)  e  7}(w,/3),  while  for  the  lateral  modes  O' =2, 4,6,8)  7}(oj,-/3)  a  -  7}(co,/3). 
These  relations  are  conveniently  summarized  by  the  following  equation 

7>,-/3)  =  (~iy+17}(6>,j3),  j  =  1,2,..., 7, 8  (21) 

Figure  3  also  shows  that  when  the  wind  direction  @w  is  not  90  degrees,  not  all  13  values  of  /3 
considered  by  SMP  are  used.  In  these  cases,  the  appearance  of  negative  values  of  /3  essentially  means 
the  double  use  of  the  values  of  Hj  and  7}  corresponding  to  that  /3.  For  example,  for  &w  =  0  degrees, 
only  values  of  Hj  and  7}  lying  between  0  and  90  degrees  are  needed.  Thus,  a  selection  and  rear¬ 
rangement  of  the  SMP  output  values  of  Hj  and  7}  are  required  for  $w  *  90  degrees. 

Once  the  proper  values  of  «  and  (3  have  been  selected,  the  wave  amplitudes  can  be  con¬ 
veniently  evaluated  by  Eq.  (9a),  where  n ,  the  wave  heading  relative  to  the  wind  direction,  is  given  in 
terms  of  /S  and  by  Eq.  (9b). 

4.2.3  Procedure  for  Arbitrary  Values  of  9 

The  Kochin  functions  Hj  are  used  as  follows  to  calculate  wave  elevations  and  slopes  at  field 
points  ( x,y )  for  cases  where  the  value  of  9,  defined  in  Eq.  (1),  does  not  correspond  exactly  to  one  of 
the  13  evenly  spaced  values  between  0  and  180  degrees  considered  by  SMP.  For  positive  values  of  0, 
i.e.,  the  field  point  has  a  positive  value  of  y,  the  value  of  Hj  is  obtained  simply  by  linear  interpolation 
between  the  two  SMP  values  of  Hj  which  bracket  the  given  value  of  0.  For  negative  values  of 
9,  Hj(-9 )  is  related  to  Hj(9)  by  noting  that  the  ship  waves  are  symmetrical  with  respect  to  9  for  the 
vertical  modes  (j- 1,3, 5, 7)  and  antisymmetric  for  the  lateral  modes,  i.e., 

Hj(-9)  =  (~iy+'Hj(9),  j  *  1,2 . 7,8  (22) 


S.  CALCULATED  WAVE  PATTERNS 

5.1  Choice  of  Ship  Hulls 

Program  SERDWAV  was  used  to  calculate  sample  wave  patterns  for  four  hulls.  Principal 
geometric  dimensions  of  the  hulls  are  given  in  Table  1.  The  table  gives  the  length  L,  beam  at  the 
free  surface  B ,  draft  H,  and  (in  the  case  of  the  twin  hull)  spacing  2Y  between  the  vertical  center- 
planes  of  the  hulls. 

Figures  4  and  5  respectively  show  the  section  contc’-s  of  the  DE1006  and  CVA59  hulls  while 
Fig.  6  shows  top  and  side  views  of  the  TAGOSTW  twin  hull.  The  DE1006  is  a  destroyer  type  hull, 
the  CVA59  represents  an  aircraft  carrier,  and  the  TAGOSTW  is  a  geometrically  simplified  version  of 
the  small-waterplane-area  twin  hull  (SWATH)  T-AGOS  described  in  [15].  The  TAGOSMH 
represents  the  monohull  version  of  the  TAGOSTW. 

5.2  Choice  of  Sea  States 

Most  of  the  runs  were  made  using  a  wind  speed  of  18  knots,  which  represents  Sea  State  4.  The 
sea  is  represented  by  five  frequencies  chosen  so  as  to  give  equal  wave  amplitudes  am,  and  13  wave 
headings  equally  spaced  90  degrees  on  either  side  of  the  wind  direction,  in  accordance  with  Eq.  (10). 
To  give  an  indication  of  the  variation  of  the  wind  speed  on  the  wave  patterns,  runs  were  also  made 
for  wind  speeds  of  12  and  42  knots,  which  correspond  respectively  to  Sea  States  2  and  8.  For  these 
three  wind  speeds.  Table  2  shows  the  five  chosen  values  of  co,  the  corresponding  wavelength  X.  and 
the  index  number  of  the  frequency  in  SMP  which  is  closest  to  the  chosen  cj.  For  the  18-knot  case. 
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Table  2  shows  that  the  shortest  wavelength  of  58  ft  is  comparable  to  the  draft  or  beam  of  most  of  the 
hulls,  while  the  longest  wavelength  of  302  ft  approximates  the  length  of  most  of  the  hulls.  For  most 
runs,  the  wind  wave  direction  was  fixed  at  45  degrees,  i.e.  quartering  seas.  Some  runs  were  also 
made  for  equal  to  0  and  90  degrees,  i.e.,  following  and  beam  seas. 

5.3  Plots  of  Wave  Patterns 

For  the  three-dimensional  and  contour  plots,  wave  elevations  and  slopes  were  calculated  over  a 
square  x  —y  grid,  with  0  x,  y  5  2 X^,  where  \M  is  the  value  of  the  largest  wavelength  X  shown 
in  Table  2  for  a  given  wind  speed.  Thus,  2  \M  is  equal  to  268,  604,  and  3,306  ft  for  wind  speeds  of 
12,  18,  and  42  knots,  respectively.  The  computation  square  was  subdivided  into  60  x  60  cells,  each 
with  a  length  equal  to  XM/30,  thus  modeling  two  cycles  of  the  maximum  wavelength.  Since  Table  2 
shows  that  the  minimum  value  of  X,  Xm,  is  approximately  equal  to  Xw/5,  this  means  that  the  shortest 
wavelength  was  modeled  by  approximately  six  points.  The  two-dimensional  plots  were  calculated  for 
the  18-knot  case  along  an  axis  x'  inclined  at  angle  a  to  the  x  axis  for  101  points  spaced  10  ft  apart, 
i.e.,  0  •£  x'  <  1000  ft. 


5.3.1  Three-Dimensional  Plots 


Figures  7a  -  7c  show  ambient  wave  elevations  for  wind  speed  V  equal  to  18  knots  (Sea  State  4) 
for  wind  direction  =  0,  45,  and  90  degrees,  respectively.  Figures  7d  -  7e  respectively  show  the 
x  and  y  wave  slopes  for  V  =  18  knots  and  =  45  degrees.  Figures  8a  -  8d  show  the  resultant 
ship  radiation  and  diffraction  waves  for  the  four  hulls  for  V  =  18  knots  and  =  45  degrees.  Fig¬ 
ures  9a  -  9d  show  the  superposition  of  the  ship  and  ambient  waves  for  the  four  hulls  for  the  same  sea 
conditions. 

Figures  7a  -  7c  shew  that  the  sea  crests  tend  to  align  perpendicular  to  the  wind  direction.  This 
is  due  to  the  fact  that  even  though  the  wind  energy  is  spread  over  a  range  of  directions,  the  principal 
energy  is  still  centered  at  the  wind  direction,  see  Eq.  (10).  Figures  7d  and  7e  show  that  the  crests  of 
the  x  and  y  wave  slopes  tend  to  align  perpendicular  to  the  x  and  y  axes,  respectively. 

Figure  8  shows  that  the  ship  waves  tend  to  have  larger  amplitudes  along  the  x-axis.  This  is 
largely  due  to  the  fact  that  the  length  of  the  ship  is  much  larger  than  its  width,  resulting  in  greater  dif¬ 
fraction  of  the  waves  traveling  along  the  length  of  the  ship.  Comparison  of  Figs.  8d  and  8e  shows 
that  the  presence  of  the  second  hull  gives  only  a  small  change  in  the  wave  pattern  compared  to  that  of 
the  monohuil. 

Comparison  of  the  plot  of  the  ambient  sea  shown  in  Fig.  7b  with  the  superposed  sea  and  ship 
waves  shown  in  Fig.  9  shows  that  for  the  three  smallest  hulls  the  presence  of  the  hull  makes  a  rela¬ 
tively  small  change  to  the  ambient  wave  pattern.  The  exception  occurs  for  the  CVA59  which  makes  a 
large  change  to  the  ambient  wave  pattern.  This  is  due  to  its  large  size  which  gives  rise  to  resultant 
large  diffraction  waves. 

5.3.2  Contour  Plots 


tion. 


A  different  view  of  the  wave  pattern  may  be  obtained  by  plotting  contours  of  equal  wave  eleva- 
Figure  10  shows  contours  of  the  ambient  sea  waves  for  the  following  cases: 


a. 

b. 

c. 

d. 

e. 


V 

V 

V 

V 

V 


-  18  knots, 
=  18  knots, 
=  18  knots, 
=  12  knots, 
=  42*knots, 


Pw 

Pw 

Pw 

Pw 

Pw 


=  0  deg 
=  45  deg 
=  90  deg 
-  45  deg 
=>  45  (leg 
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Figures  1  la  -  d  show  the  contour  lines  of  the  superposed  ship  and  ambient  waves  for  the  four  hulls 
for  the  wind  conditions  given  in  Fig.  10b.  Figures  lie  and  f  respectively  show  the  superposed  com 
tour  lines  for  the  DE1006  for  the  wind  conditions  given  in  Figs.  lOd  and  e. 

As  was  previously  noticed  in  the  case  of  the  three-dimensional  plots,  Fig.  10  clearly  shows  that 
the  contour  lines  tend  to  align  perpendicular  to  the  wind  direction.  In  the  absence  of  any  energy 
spread,  the  contours  would,  of  course,  be  straight  lines  exactly  perpendicular  to  the  wind  direction. 

Figures  10b,  d,  and  e  show  the  interesting  feature  that  the  contour  lines  for  a  given  value  of 
45  deg),  the  ambient  wave  contour  lines  for  different  wind  speeds  are  quite  similar.  This  is 
largely  a  consequence  of  the  size  of  the  computation  grid  scaled  to  be  equal  to  twice  the  maximum 
wavelength,  as  well  as  the  similar  ratios  between  the  five  component  wavelengths  which  are  selected 
to  model  the  sea  (see  Table  2). 

Figures  lla-d  again  show  that  for  the  three  smallest  hulls,  the  presence  of  the  hull  in  a  18-knot 
sea  makes  a  relatively  small  change  to  the  contour  lines  of  the  ambient  waves  alone,  shown  in  Fig. 
10b.  On  the  other  hand,  as  previously  noted  in  the  three-dimensional  plots,  the  presence  of  the  990-ft 
CVA59  makes  a  large  change  to  the  contour  lines.  This  effect  of  ship  size  relative  to  ambient 
wavelengths  on  the  resultant  perturbation  to  the  ambient  waves  is  also  clearly  illustrated  in  Figs,  lie 
and  f.  Figure  lie  shows  that  the  DEI 006  hull  makes  a  large  change  to  the  contours  of  the  short- 
wavelength  12-knot  sea  shown  in  Fig.  lOd.  On  the  other  hand,  Fig.  Ilf  shows  that  this  hull  hardly 
makes  any  changes  to  the  contours  of  the  long-wavelength  42-knot  sea,  shown  in  Fig.  lOe. 

5.3.3  Two-Dimensional  Plots 

Figures  12  a-c  show  curves  of  the  ambient  waves  and  superposed  ship  and  ambient  waves  for 
the  DE1QC6  in  the  presence  of  a  V  -  18  knot,  =  0  deg  ambient  sea  along  an  axis  x'  inclined 
respectively  at  0,  45,  and  90  degrees  to  the  x  axis.  Figures  12d-f  and  12g-i  show  corresponding 
results  for  @w  —  45  and  90  deg,  respectively.  Figures  12j-l  show  curves  of  the  ambient  and  super¬ 
posed  waves  along  the  x  axis  (a  =  0  deg)  for  the  CVA59  for  V  -  18  knots  and  0W  =  0,  45,  and  90 
deg,  respectively. 

As  previously  noted  in  the  three-dimensional  plots.  Fig.  12  shows  that  the  largest  ship  waves 
occur  along  the  x  axis,  a  =  0  deg,  due  to  the  large  diffraction  of  waves  in  this  direction.  For  the 
DE1006,  Figs.  12a,  d,  and  g  show  that  the  ship  waves  are  relatively  small  but  still  noticeable  at  x  = 
1000  ft.  The  figure  shows  that  the  ship  waves  for  a  =  45  and  90  deg  are  smaller  than  those  for 
a  =  0  deg.  Figures  12c,  f,  and  i  show  the  expected  trend  that  the  ship  waves  for  a  =  90  deg 
increase  as  the  wind  direction  shifts  from  0  to  90  deg.  Thus,  Figs.  12g-i  show  that  for  =  90  deg, 
the  ship  waves  for  all  values  of  a  are  still  noticeable  (relative  to  the  ambient  waves)  at  x  '  =  1000  ft. 

Figures  12j-I  show  the  much  larger  magnitude  of  the  ship  waves  in  the  case  of  the  CVA59.  For 
example,  at  x  —  1000  ft,  the  ship  waves  are  still  larger  than  the  ambient  waves. 

6.  SUMMARY 

This  report  has  discussed  in  detail  the  procedure  for  generating  spatial  realizations  of  the  far 
field  ship  radiation  and  diffraction  waves  in  the  presence  of  random  ambient  sea  waves.  The  pro¬ 
cedure  essentially  consists  of  generalizing  a  previously  described  technique  for  obtaining  the  ship 
wave  pattern  due  to  a  single  incident  wave,  using  the  Kochin  function  approach,  to  the  multiple  waves 
of  a  random  sea.  The  random  sea  waves  are  expressed  as  a  double  sum  over  frequency  and  heading 
of  randomly  phased  sinusoidal  components.  The  ship  waves  are  expressed  as  a  triple  sum  over  eight 
ship  wave  modes  and  the. above  frequencies  and  headings.  The  eight  wave  modes  consist  of  radiation 
wave  modes  due  to  the  six  oscillatory  motions  of  the  ship,  and  the  vertical  and  lateral  components  of 
the  diffraction  wave  due  to  scattering  by  the  hull  surface.  It  is  shown  that  by  a  combination  of  inter¬ 
polation,  and  use  of  the  symmetry  and  antisymmetry  properties  of  the  ship  wave  modes,  the  discrete 
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set  of  Kochin  functions  and  motion  transfer  functions  output  by  the  Ship  Motion  Program  (SMP)  may 
be  used  to  calculate  ship  wave  elevations  for  random  sea  waves  at  arbitrary  field  point  locations. 

It  is  shown  that  the  horizontal  x  and  y  slopes  of  the  waves  may  be  conveniently  obtained  from 
the  expression  for  wave  elevation  by  multiplication  by  the  wavenumbers  in  the  x  and  y  directions, 
respectively.  The  formulation  for  a  monohull  is  also  generalized  to  the  twin  hull  case  by  considering 
the  vertical  plane  midway  between  the  hulls  as  a  plane  of  symmetry.  It  is  shown  that  this  approach 
may  qualitatively  give  the  wave  pattern  for  field  points  outside  the  channel  area  between  the  two 
hulls. 


Sea  and  ship  waves  are  presented  for  four  ship  hulls  in  the  presence  of  various  ambient  seas  in 
the  form  of  three-dimensional,  two-dimensional,  and  contour  plots.  Three  of  the  hulls  are  250  to  308 
feet  long,  while  the  fourth  hull  is  the  990-ft  aircraft  carrier  CVA59.  The  plots  for  the  sea  waves 
show  that  their  contours  tend  to  be  aligned  perpendicular  to  the  wind  direction.  The  ship  radiation 
and  diffraction  waves  are  largest  along  x ,  the  direction  parallel  to  the  ship  longitudinal  axis.  The  ship 
waves  along  y ,  the  horizontal  direction  perpendicular  to  x ,  increase  as  the  wind  direction  shifts  from 
x  to  y .  In  a  Sea  State  4,  the  presence  of  the  250  and  308  ft  ships  makes  a  moderately  small  change 
to  the  ambient  waves,  while  the  990-ft  CVA59  makes  a  large  change.  Similarly,  relative  to  the  sea 
waves,  the  ship  waves  due  to  the  308-ft  DE1006  are  large  in  a  Sea  State  2  and  nearly  imperceptible  in 
a  Sea  State  8.  The  presence  of  a  second  hull  gives  a  small,  though  perceptible,  change  to  the  wave 
pattern  of  a  monohull. 
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Table  1  —  Hull  Geometric  Characteristics 


Ship  Hull 

L(ft) 

Bfft) 

H(ft) 

2Y(ft) 

1.  DEI 006 

□ 

35.9 

12.1 

— 

2.  CVA59 

1 

129.3 

35.8 

— 

3.  TAGOSTW 

250. 

7.0 

24.3 

80. 

4.  TAGOSMH 

250. 

7.0 

24.3 

— 

Table  2  —  Wave  Components 
for  Sea  States  2,  4,  and  8 


o.'(rad/sec) 

X(ft) 

SMP# 

Sea  State  2 

,  V  =  12 

knots 

1.229 

134 

23 

1.433 

99 

25 

1.635 

76 

26 

1.922 

55 

28 

2.774 

26 

30 

Sea  State  4 

< 

11 

00 

knots 

0.818 

302 

17 

0.956 

221 

20 

1.092 

170 

22 

1.285 

123 

24 

1.874 

58 

27 

Sea  State  8 

II 

> 

knots 

0.349 

1,653 

4 

0.410 

1,204 

5 

0.468 

924 

8 

0.551 

666 

11 

0.806 

312 

17 
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—  Relation  of  sea  wave  headings  to  wind  direction 


Section  contours  for  carrier  hull  form  CVA59 


SIDE  VIEW  OF  TAGOSTW  TWIN  HULL 

Figure  6  —  Geometry  of  small  waterplane  twin  hull  TAGOSTW 


(a)  Elevation,  V=  18  knots,  !iw  =  0  deg 


(b)  Elevation,  V  =  18  knots,  —  45  deg 
Figure  7  —  Three-dimensional  plots  of  ambient  wave  elevations  and  slopes 
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(c)  Elevation,  V=18  knots,  (jW  =  90  deg 


(d)  x-slope,  V=  18  knots,  (3W  =  45  deg 
Fig.  7  (Continued)  —  Three-dimensional  plots  of  ambient  wave  elevations  and  slopes 
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(e)  y-slope,  V=18  knots,  {3W  —  45  deg 


Fig.  7  (Continued)  —  Three-dimensional  plots  of  ambient  wave  elevations  and  slopes 
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(a)  Destroyer  hull  form  DEI 006 


(b)  Carrier  hull  form  CVA59 

Figure  8  —  Ship  radiation  and  diffraction  waves  for  V=18  knots,  =  45  deg 
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(c)  Small  waterplane  area  twin  hull  TAGOSTW 


(d)  Small  waterplane  area  monohul!  TAGOSMH 
Fig.  8  (Continued)  —  Ship  radiation  and  diffraction  waves  for  18  knots,  (3W  ■  45  deg 
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(a)  Destroyer  hull  form  DEI 006 


(b)  Carrier  hull  form  CVA59 

Figure  9  —  Superposed  ship  and  ambient  waves  for  V=18  knots,  (3W  ~  45  deg 
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(c)  Small  waterplane  area  twin  hull  TAGOSTW 


(d)  Small  waterplane  area  monohull  TAGOSMH 
Fig.  9  (Continued,  —  Superposed  ship  and  ambient  waves  for  V=  18  knots,  -  45  deg 
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(b)  V  =  18  knots,  0  w  -  45  deg 


Figure  10  —  Contour  plots  of  ambient  wave  elevations 
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(d)  V  =  12  knots,  =  45  deg 


Fig.  10  (Continued)  —  Contour  plots  of  ambient  wave  elevations 
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(e)  V  =42  knots,  fiw  =  45  deg 


Fig.  10  (Continued)  —  Contour  plots  of  ambient  wave  elevations 


(a)  DE1006,  V=18  knots,  0W  =  45  deg 


1 1  —  Contour  plots  of  superposed  ship  and  ambient  wave  elevations 


(d)  TAGOSMH,  V-18  knots,  0W  »  45  deg 
Fig.  1 1  (Continued)  —  Contour  plots  of  superposed  ship  and  ambient  wave  elevations 
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(e)  DE1006,  V= 12  knots,  -  45  deg 


(0  DE1006,  V =42  knots,  =  45  deg 


Fig.  1 1  (Continued)  —  Contour  plots  of  superposed  ship  and  ambient  wave  elevations 
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WAVE  ELEVATION  (FT)  WAVE  ELEVATION  (FT) 


AXIAL  DISTANCE  (FT) 

(b)  DE1006,  a  =  45  deg,  V  =  18  knots,  0  w  =  0  deg 

Figure  12  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
'and  ambient  wave  elevations  along  fixed  directions 
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WAVE  ELEVATION  (FT) 


(c)  DE1006,  a  =  90  deg,  V=  18  knots,  0  w  =  0  deg 


Fig.  12  (Continued)  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
and  ambient  wave  elevations  along  fixed  directions 


WAVE  ELEVATION  (FT)  WAVE  ELEVATION  (FT) 


AXIAL  DISTANCE  (FT) 


(e)  DE1006,  a  -  45  deg,  V=18  knots,  0  w  -  45  deg 


0.0  200.0  400.0  600.0  800.0  1000.0 

AXIAL  DISTANCE  (FT) 

(f)  DE1006,  a  =  90  deg,  V=  18  knots,  0  w  =45  deg 

Fig.  12  (Continued)  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
and  ambient  wave  elevations  along  fixed  directions 


WAVE  ELEVATION  (FT)  WAVE  ELEVATION  (FT) 


(g)  DE1006,  a  -  0  deg,  V  =  18  knots,  /3  w  -  90  deg 


AXIAL  DISTANCE  (FT) 

(h)  DE1006,  a  -  45  deg,  V  =  18  knots,  (3  w  =  90  deg 


Fig.  12  (Continued)  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
and  ambient  wave  elevations  along  fixed  directions 
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WAVE  ELEVATION  (FT)  WAVE  ELEVATION  (FT) 


AXIAL  DISTANCE  (FT) 


(i)  DE1006,  a  =  90  deg,  V=18  knots,  0  w  =  90  deg 


(j)  CVA59,  a  -  0  deg,  V=  18  knots,  0  &  =  0  deg 
•  • 

Fig.  12  (Continued)  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
and  ambient  wave  elevations  along  fixed  directions 
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WAVE  ELEVATION  (FT)  WAVE  ELEVATION  (FT) 


(k)  CVA59,  a  =  0  deg,  V=  18  knots,  0  w  =  45  deg 


AXIAL  DISTANCE  (FT) 

(1)  CVA59,  a  -  0  deg,  V=  18  knots,  0  w  »  90  deg 

Fig.  12  (Continued)  —  Two-dimensional  plots  of  ambient,  and  superposed  ship 
and  ambient  wave  elevations  along  fixed  directions 
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